Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution          下载免费PDF全文

Wencheng Jin  Hao Xu  Chloé Arson  Seth Busetti 《国际地质力学数值与分析法杂志》2017,41(2):223-250

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

The influence of coupled physical swelling and chemical reactions on deformable geomaterials     

Yue Ma  Xiao‐Hui Chen  Lee J. Hosking  Hai‐Sui Yu  Hywel R. Thomas  Simon Norris 《国际地质力学数值与分析法杂志》2021,45(1):64-82

Coupled thermo‐hydro‐mechanical‐chemical modelling has attracted attention in past decades due to many contemporary geotechnical engineering applications (e.g., waste disposal, carbon capture and storage). However, molecular‐scale interactions within geomaterials (e.g., swelling and dissolution/precipitation) have a significant influence on the mechanical behaviour, yet are rarely incorporated into existing Thermal‐Hydro‐Mechanical‐Chemical (THMC) frameworks. This paper presents a new coupled hydro‐mechanical‐chemical constitutive model to bridge molecular‐scale interactions with macro‐physical deformation by combining the swelling and dissolution/precipitation through an extension of the new mixture‐coupling theory. Entropy analysis of the geomaterial system provides dissipation energy, and Helmholtz free energy gives the relationship between solids and fluids. Numerical simulation is used to compare with the selected recognized models, which demonstrates that the swelling and dissolution/precipitation processes may have a significant influence on the mechanical deformation of the geomaterials.  相似文献   

Quantifying flows along hydrological pathways by applying a new filtering algorithm in conjunction with master recession curve analysis          下载免费PDF全文

Ronan J. O'Brien  Bruce D. Misstear  Laurence W. Gill  Paul M. Johnston  Raymond Flynn 《水文研究》2014,28(26):6211-6221

Quantifying the proportion of the river hydrograph derived from the different hydrological pathways is essential for understanding the behaviour of a catchment. This paper describes a new approach using the output from master recession curve analysis to inform a new algorithm based on the Lyne and Hollick ‘one‐parameter’ signal analysis filtering algorithm. This approach was applied to six catchments (including two subcatchments of these) in Ireland. The conceptual model for each catchment consists of four main flow pathways: overland flow, interflow, shallow groundwater and deep groundwater. The results were compared with those of the master recession curve analysis, a recharge coefficient approach developed in Ireland and the semi‐distributed, lumped and deterministic hydrological model Nedbør‐Afstrømings‐Model. The new algorithm removes the ‘free variable’ aspect that is typically associated with filtering algorithms and provides a means of estimating the contribution of each pathway that is consistent with the results of hydrograph separation in catchments that are dominated by quick response pathways. These types of catchments are underlain by poorly productive aquifers that are not capable of providing large baseflows in the river. Such aquifers underlie over 73% of Ireland, ensuring that this new algorithm is applicable in the majority of catchments in Ireland and potentially in those catchments internationally that are strongly influenced by the quick‐responding hydrological pathways. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Is there a superior conceptual groundwater model structure for baseflow simulation?          下载免费PDF全文

Michael Stoelzle  Markus Weiler  Kerstin Stahl  Andreas Morhard  Tobias Schuetz 《水文研究》2015,29(6):1301-1313

Previous work has shown that streamflow response during baseflow conditions is a function of storage, but also that this functional relationship varies among seasons and catchments. Traditionally, hydrological models incorporate conceptual groundwater models consisting of linear or non‐linear storage–outflow functions. Identification of the right model structure and model parameterization however is challenging. The aim of this paper is to systematically test different model structures in a set of catchments where different aquifer types govern baseflow generation processes. Nine different two‐parameter conceptual groundwater models are applied with multi‐objective calibration to transform two different groundwater recharge series derived from a soil‐atmosphere‐vegetation transfer model into baseflow separated from streamflow data. The relative performance differences of the model structures allow to systematically improve the understanding of baseflow generation processes and to identify most appropriate model structures for different aquifer types. We found more versatile and more aquifer‐specific optimal model structures and elucidate the role of interflow, flow paths, recharge regimes and partially contributing storages. Aquifer‐specific recommendations of storage models were found for fractured and karstic aquifers, whereas large storage capacities blur the identification of superior model structures for complex and porous aquifers. A model performance matrix is presented, which highlights the joint effects of different recharge inputs, calibration criteria, model structures and aquifer types. The matrix is a guidance to improve groundwater model structures towards their representation of the dominant baseflow generation processes of specific aquifer types. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Including spatial distribution in a data‐driven rainfall‐runoff model to improve reservoir inflow forecasting in Taiwan          下载免费PDF全文

Meng‐Jung Tsai  Robert J. Abrahart  Nick J. Mount  Fi‐John Chang 《水文研究》2014,28(3):1055-1070

Multi‐step ahead inflow forecasting has a critical role to play in reservoir operation and management in Taiwan during typhoons as statutory legislation requires a minimum of 3‐h warning to be issued before any reservoir releases are made. However, the complex spatial and temporal heterogeneity of typhoon rainfall, coupled with a remote and mountainous physiographic context, makes the development of real‐time rainfall‐runoff models that can accurately predict reservoir inflow several hours ahead of time challenging. Consequently, there is an urgent, operational requirement for models that can enhance reservoir inflow prediction at forecast horizons of more than 3 h. In this paper, we develop a novel semi‐distributed, data‐driven, rainfall‐runoff model for the Shihmen catchment, north Taiwan. A suite of Adaptive Network‐based Fuzzy Inference System solutions is created using various combinations of autoregressive, spatially lumped radar and point‐based rain gauge predictors. Different levels of spatially aggregated radar‐derived rainfall data are used to generate 4, 8 and 12 sub‐catchment input drivers. In general, the semi‐distributed radar rainfall models outperform their less complex counterparts in predictions of reservoir inflow at lead times greater than 3 h. Performance is found to be optimal when spatial aggregation is restricted to four sub‐catchments, with up to 30% improvements in the performance over lumped and point‐based models being evident at 5‐h lead times. The potential benefits of applying semi‐distributed, data‐driven models in reservoir inflow modelling specifically, and hydrological modelling more generally, are thus demonstrated. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Quantifying the impact of landscape changes on hydrological variables in the alpine and cold region using hydrological model and remote sensing data     

Zizhen Jin  Qiudong Zhao  Xiang Qin  Jingtian Zhang  Hui Zhang  Jia Qin  Yu Qin  Hongyuan Li  Jizu Chen  Yushuo Liu  Yanzhao Li  Lihui Wang 《水文研究》2021,35(10):e14392

Quantifying the impact of landscape on hydrological variables is essential for the sustainable development of water resources. Understanding how landscape changes influence hydrological variables will greatly enhance the understanding of hydrological processes. Important vegetation parameters are considered in this study by using remote sensing data and VIC-CAS model to analyse the impact of landscape changes on hydrology in upper reaches of the Shule River Basin (URSLB). The results show there are differences in the runoff generation of landscape both in space and time. With increasing altitude, the runoff yields increased, with approximately 79.9% of the total runoff generated in the high mountains (4200–5900 m), and mainly consumed in the mid-low mountain region. Glacier landscape produced the largest runoff yields (24.9% of the total runoff), followed by low-coverage grassland (LG; 22.5%), alpine cold desert (AL; 19.6%), mid-coverage grassland (MG; 15.6%), bare land (12.5%), high-coverage grassland (HG; 4.5%) and shrubbery (0.4%). The relative capacity of runoff generation by landscapes, from high to low, was the glaciers, AL, LG, HG, MG, shrubbery and bare land. Furthermore, changes in landscapes cause hydrological variables changes, including evapotranspiration, runoff and baseflow. The study revealed that HG, MG, and bare land have a positive impact on evapotranspiration and a negative impact on runoff and baseflow, whereas AL and LG have a positive impact on runoff and baseflow and a negative impact on evapotranspiration. In contrast, glaciers have a positive impact on runoff. After the simulation in four vegetation scenarios, we concluded that the runoff regulation ability of grassland is greater than that of bare land. The grassland landscape is essential since it reduced the flood peak and conserved the soil and water.  相似文献   

珠江三角洲地区登革热流行风险空间模拟与预测     

郑斓  任红艳  施润和  鲁亮 《地球信息科学学报》2019,21(3):407-416

近年来日益严重的登革热疫情已在中国南部地区形成疫情高发区,并对中国的公共卫生安全形成了一定的威胁。登革热主要受到区域内复杂的自然环境条件以及社会经济因素的影响,而利用地理空间分析方法和模型探究登革热疫情的影响因素,并对其未来流行风险的空间分布进行模拟,是有效开展登革热预防控制工作的重要基础。本文收集了珠江三角洲地区2010-2014年的登革热病例资料和土地利用、人口密度两种社会经济要素数据,构建土地利用回归（LUR）模型以分析登革热疫情与不同空间范围内的土地利用和人口密度之间的关系,并结合SLEUTH模型获取的2030年土地利用数据以及基于人口密度预测模型获取的2030年人口密度数据,预测珠江三角洲地区2030年登革热疫情风险的空间分布。结果表明,社会经济要素对登革热疫情空间分布的影响在不同范围内存在差异,半径分别为10、7、10、2和1 km的缓冲区内的人口密度、草地、城镇用地、林地和耕地进入LUR模型并对疫情有显著的影响（相关系数分别为0.779、-0.473、0.818、-0.642和-0.403）,所构建的LUR模型效果较好（调整R²为0.796,F=390.409,P<0.01）,留一交叉检验结果显示模型的相对均方根误差为0.7046,预测值与实测值的拟合精度达到0.7101。2030年城市空间扩展的区域主要分布在深圳、东莞以及广佛的交界地区,而登革热风险预测模型表明2030年登革热疫情风险较大的区域与珠江三角洲城镇用地占比、人口分布较高的地区有高度的一致性,尤其是广佛地区。因此,LUR模型可以较好地预测登革热疫情的空间分布,从而为当地卫生部门防控登革热提供方法支持。  相似文献   

Ekman模型时间变化风应力系数的伴随参数估计     

李国庆  高艳秋  张继才 《海洋与湖沼》2019,50(5):979-993

本文基于时间分布参数设置,利用伴随同化方法,反演了Ekman模型中随时间变化的风应力拖曳系数,并在孪生实验和实际实验中对该方法进行了验证。在孪生实验中,研究了参数反演结果对不同影响因素的响应,包括:风速分布、风应力系数分布、风应力系数初始猜测值、风应力系数独立变量个数、观测数据误差和观测的深度。孪生实验结果验证了伴随同化方法反演Ekman模型中时变风应力系数的有效性,具体包括如下五个方面结论:1)不同风速分布下均能成功反演出不同风应力拖曳系数分布; 2)反演结果对初始猜测值较为敏感,风应力系数初始猜测值越接近给定值,反演结果越好;3)风应力系数独立点个数的选取会显著影响反演结果,合理的选择有利于提高反演效率及减小观测数据误差;4)观测误差能够影响反演结果,观测数据误差在20%以下时能取得合理的反演结果; 5)反演结果对观测数据的表层和次表层流速更为敏感,这是由Ekman流的物理性质决定的。实际实验,利用百慕大锚系试验平台的风速和流速数据,去除周期性潮流和地转流成分后得到Ekman流成分,并作为观测输入到该同化模型,反演出了适用于该区域和该时段的随时间变化的风应力系数。通过比较模拟流速和观测流速,证明利用伴随同化方法能从实测数据中反演出合理的时变风应力系数,对于海洋模型风应力系数的确定是一项有益的尝试。  相似文献   

Severe soil erosion during a 3‐day exceptional rainfall event: combining modelling and field data for a fallow cereal field          下载免费PDF全文

Manuel López‐Vicente  Laura Quijano  Leticia Gaspar  Leticia Palazón  Ana Navas 《水文研究》2015,29(10):2358-2372

Exceptional rainfall events cause significant losses of soil, although few studies have addressed the validation of model predictions at field scale during severe erosive episodes. In this study, we evaluate the predictive ability of the enhanced Soil Erosion and Redistribution Tool (SERT‐2014) model for mapping and quantifying soil erosion during the exceptional rainfall event (~235 mm) that affected the Central Spanish Pyrenees in October 2012. The capacity of the simulation model is evaluated in a fallow cereal field (1.9 ha) at a high spatial scale (1 × 1 m). Validation was performed with field‐quantified rates of soil loss in the rills and ephemeral gullies and also with a detailed map of soil redistribution. The SERT‐2014 model was run for the six rainfall sub‐events that made up the exceptional event, simulating the different hydrological responses of soils with maximum runoff depths ranging between 40 and 1017 mm. Predicted average and maximum soil erosion was 11 and 117 Mg ha^?1 event^?1, respectively. Total soil loss and sediment yield to the La Reina gully amounted to 16.3 and 9.0 Mg event^?1. These rates are in agreement with field estimations of soil loss of 20.0 Mg event^?1. Most soil loss (86%) occurred during the first sub‐event. Although soil accumulation was overestimated in the first sub‐event because of the large amount of detached soil, the enhanced SERT‐2014 model successfully predicted the different spatial patterns and values of soil redistribution for each sub‐event. Further research should focus on stream transport capacity. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Modified Rindler acceleration as a nonlinear electromagnetic effect     

《Astroparticle Physics》2015

The model proposed originally by Mannheim and Kazanas for fitting the shapes of galactic rotation curves has recently been considered by Grumiller to describe gravity of a central object at large distances. Herein we employ the same geometry within the context of nonlinear electrodynamics (NED). Pure electrical NED model is shown to generate the novel Rindler acceleration term in the metric which explains anomalous behaviors of test particles/satellites. Remarkably a pure magnetic model of NED yields flat rotation curves that may account for the missing dark matter. Weak and strong energy conditions are satisfied in such models of NED.  相似文献